Urban traffic from the perspective of dual graph 



Mao-Bin Hu 1 - 2 Q Rui Jiang 1 , Yong-Hong Wu 2 , Wen-Xu Wang 3 , and Qing -Song Wufl 
1 School of Engineering Science, University of Science and Technology of China, Hefei 230026, P.R.China 
2 Department of Mathematics and Statistics, Curtin University of Technology, Perth WA6845, Australia 
3 Department of Electronic Engineering, City University of Hong Kong, Hong Kong SAR, P.R. China 

(Dated: February 2, 2008) 

In this paper, urban traffic is modeled using dual graph representation of urban transportation 
network where roads are mapped to nodes and intersections are mapped to links. The proposed 
model considers both the navigation of vehicles on the network and the motion of vehicles along 
roads. The road's capacity and the vehicle-turning ability at intersections are naturally incorporated 
in the model. The overall capacity of the system can be quantified by a phase transition from free 
flow to congestion. Simulation results show that the system's capacity depends greatly on the 
topology of transportation networks. In general, a well-planned grid can hold more vehicles and its 
overall capacity is much larger than that of a growing scale-free network. 

PACS numbers: 45.70.Vn, 89.75.Hc, 05.70.Fh 



I. INTRODUCTION 

Traffic problem is of great importance for the safety 
and convenience of modern society [l|, 0, H, 0, H, [if. 
The traffic research was mainly focused on the highway 
traffic El 0, El or the traffic on a well-planned grid lat- 
tice [HTH, [g]. Recently, empirical evidences have shown 
that many transportation systems can be described by 
complex networks characterized by the small world |7[ 
and/or scale- free properties M. The prototypes include 
urban road network [9, [l(| [U 0, EH > nation-wide road 
network [T3, public transport network [TBI. HHl, railway 
network [171 . Il8j , airline transport systems [19(, and so 
on. And due to the importance of large communication 
networks such as the Internet and WWW, much inves- 
tigation have been focused on ensuring free traffic flow 
and avoiding traffic congestion on networks. 

In urban traffic, one natural approach is the "primal" 
representation which considers intersections as nodes and 
road segments as links between nodes. However, this 
representation can not characterize the common notion 
of "main roads" in urban traffic systems. Moreover, it 
violates the intuitive notion that an intersection is where 
two roads cross, not where four roads begin, so that it 
can not represent the way we usually tell each other how 
to navigate the road network [l4j], e.g., "stay in Manning 
Road ignoring 4 crossings, until you reach the Kent Street 
and turn right." 

Here, we study the urban traffic by employing the dual 
representation of the road network, in which a node rep- 
resents a single road, and two nodes are linked if their 
corresponding roads intersect. This kind of transforma- 
tion was firstly proposed in the field of urban pla nning 
and design with the name Space Syntax [2(| HH, [22j , and 
has been used recently to study the topological properties 
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of urban roads H, E3, ED, El El- In this dual presenta- 
tion, the "main roads" of the city can be represented by 
the "hub nodes" in the network. The node degree is not 
limited and it was found that the dual degree distribution 
of some cities follows a power-law. 

One shortcoming of the dual representation is the 
abandonment of metric distance so that a street is one 
node no matter how long it is. However, the metric dis- 
tance is the core of most traffic flow studies. To avoid this 
problem, this paper propose a traffic flow model consid- 
ering the movement of vehicles along the road. We show 
that by using two parameters characterizing the road and 
intersection capacity, it is possible to simulate the ur- 
ban traffic and seek out the overall capacity of the urban 
transportation system. 



II. DUAL GRAPH OF URBAN NETWORK 
AND THE TRAFFIC MODEL 

The basic steps of dual representation can be explained 
as following: one should firstly define which road seg- 
ments belong together. Previous studies grou ped the seg- 
ments by line of sight by a driver [2(| l2ll |22|] , their street 
name [9|, L!J| , or by using a threshold on the angle of in- 
cidence of segments at an intersection 

[Hill, El- Then, 

in the derived "dual graph" or the "connectivity graph" , 
each road is turned into one node, while each intersection 
is turned into one link. 



A. Characterize Dual Graph Considering Traffic 
Problem 

With this new paradigm, one can look at urban traffic 
problem from a new perspective. From common sense, 
the "main roads" can usually be characterized by long 
length, wide road width, many intersections, and high 
traffic efficiency. When the main roads are congested, the 
entire transportation system will be in danger of capacity 
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drop. While the system will remain high efficiency even 
when some minor roads are jam. 

In dual network, the major roads in urban system can 
be represented with hub nodes with many links, large 
capacity, and more handling ability. In particular, when 
considering the traffic flow problem, the characteristics 
of nodes can be introduced in the dual graph of urban 
network: 

(1) Degree fc, and degree distribution: The degree (or 
connectivity) of node i is the number of links connecting 
to that node. In the dual graph of urban network, the 
degree corresponds to the number of intersections along 
road i. Thus there is no limit of node degree and so the 
"main roads" can be represented by the "hub nodes" in 
the system. 

The way the degree is distributed among the nodes 
can be investigated by calculating the degree distribution 
P(k), i.e., the probability of finding nodes with k links. 
Networks with a power-law distribution are called scale- 
free [|. 

(2) Capacity C{\ The capacity of road i is the maxi- 
mum number of vehicles the road can hold. This value 
can be calculated by: Ci = Li/l v , where Li is the length 
of the road, and l v is the vehicle's average length. In 
the simulation, we assume that the length of each road 
segment are same, thus Ci is proportional to the degree 
of the node: Ci — a x fcj, where a denotes the maximum 
number of vehicles that one road segment can hold. The 
system's total capacity is the sum of the capacity of all 
node, that is C t — T,id. 

(3) Turning Ability Tf. The turning ability is the max- 
imum number of vehicles turning from the road to neigh- 
boring roads per time step. This value reflects the ca- 
pacity of intersections along the road, and can be also 
measured in real practice. In the dual representation 
of urban traffic, Ti is the maximum number of vehicles a 
node can send to neighboring nodes. Without losing gen- 
erality, we assumed that all intersections can handle the 
same number of vehicle-turning from each road, thus Ti 
is also proportional to the degree of the node: T = fixki, 
where (3 denotes the ability of one intersection. 



B. Traffic Model 

In dual network, the trajectory of a vehicle can be 
interpreted as traveling along some roads (nodes) for 
some distance, and jumping from a node to another node 
through an link representing an intersection. The pro- 
posed traffic model considers both the navigation of ve- 
hicles in the network and the movement of vehicles along 
roads. The model ispartially motivated by the vehicular 
traffic models [l], [U, Hj and by the traffic models of packet 
flow on the Internet 0, 0, [H, H3, HI • On the base of 
underlying dual infrastructure (discussed in Section III 
and IV), the system evolves in parallel according to the 
following rules: 

1. Add Vehicles - At each time step, vehicles are added 



to the system with a given rate R at randomly selected 
nodes and each new vehicle is given a random destination 
node. 

2. Navigate Vehicles - If a vehicle's destination is found 
in its nearest neighborhood, its direction will be set to the 
target. Otherwise, its direction will be set to a neighbor 
n with probability: 

L<t> 

where the sum runs over all neighbors, and 4> is an ad- 
justable parameter reflecting the driver's preference with 
the roads. It is assumed that the vehicles use a local 
routing strategy: they are more ready to use a neighbor- 
ing "main road" when <j) > 0, and they are more likely to 
go to a minor road when <fi < 0. Once a vehicle reach its 
destination, it will be removed from the system. 

3. Vehicles Motion along Roads - The intersections 
along a road i are numbered with serial integers from 1 
to ki. We use these integers to reflect the sequence of 
intersections along the road. When a vehicle enters a 
road at the mth intersection and leaves at the nth inter- 
section, it has to travel the distance of d = 1$ x \m — n\ 
along the road, where lo denotes the length of one road 
segment, and |...| denotes taking the absolute value. For 
the velocity, we assumed that the traffic flow is homoge- 
neous along one given road. In each time step, the mean 
velocity of vehicles on the road is calculated using the 
following equation: 

Vi = V max x (1.0 - p^, (2) 

where V max is the maximum velocity in the urban system, 
and pi = Ni/ Li is the local vehicle density on that road. 
Equation @ simply approximates that the vehicles will 
take V m ax when there is no vehicle on the road, and the 
velocity will decreases linearly until no vehicle can move 
when p = 1.0. In each time step, the distance of each 
vehicle in the road will decrease by the amount of Vi until 
d = 0, i.e., it reaches the preferred exiting intersection. 

4. Vehicle- Turning at Intersections - At each step, only 
at most Ti vehicles can be sent from a node for the neigh- 
boring nodes. When the number of the vehicles at a se- 
lected node is full (see the introduction of Ci), the node 
won't accept any more vehicles and the vehicle will wait 
for the next opportunity. 

One can see that the core parameters characterizing 
the transportation system are a, (3 and V max . And pa- 
rameter <f> characterizes the behavior of drivers. In the 
simulations, without losing generality, the road segment 
length Iq is assumed to be 500 meters, and the average 
vehicle length l v is 10 meters. Thus we set a = 50, 
i.e., one road segment between two successive intersec- 
tions can hold 50 vehicles at most. The parameter j3 
is also set to be constant for all intersections, meaning 
that the vehicle-turning ability for each intersection in 
the system is the same. And each time step is assumed 
to represent 5 seconds in reality. Thus if the maximum 
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FIG. 3: (Color online)J? c vs Vmax with different values of f3 . 
The data are obtained by averaging R c with ten simulations. 



FIG. 1: (Color online) Typical evolution of vehicle number 
in the grid system. The critical generating rate of vehicles is 
R c = 113. 
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FIG. 2: (Color online) R c vs /3 with different values of Vmax- 
The data are obtained by averaging R c with ten simulations. 

velocity is 20m/s, the maximum velocity in Eql^will be 
Vm ax = 100. With changing the values of (3 and V max , 
different ability of intersections and different road condi- 
tions can be simulated. 



III. SIMULATION OF URBAN TRAFFIC ON A 
WELL-PLANNED LATTICE NETWORK 

In this section, the traffic on a well-planned grid ur- 
ban network is simulated. The transportation network 
consists of 24 north-south roads crossing with 24 east- 
west roads. The studied system can be represented in 
dual graph by 48 nodes each with the same connectivity 



k = 24. Thus there is no topological "main roads" in this 
system. For the case, 4> — is applied for all vehicles, i.e., 
they randomly select direction from neighboring roads if 
its destination is not found. Figure Q] displays the typical 
evolution of N c (t), i.e., the number of vehicles within the 
system for different generating rate R. One can see that 
when R < R Cl N c will increase first and then come to sat- 
uration, indicating the balance of the number of vehicles 
entering the system and the number of vehicles reaching 
their destinations. However, when R > R c , N c can not 
remain constant. It will suddenly increase and quickly 
reach the the system's total capacity. Thus the system 
is congested and the vehicles accumulate in the system. 
Therefore, the critical generating rate R c can be used to 
characterize the phase transition from free flow state to 
congestion. And the system's overall traffic capacity can 
be measured by R c above which the system will enter the 
congested state. 

To better understand the effects of intersection capac- 
ity and road condition on the overall efficiency, simula- 
tions with different values of [3 and V ma x are carried out. 
Figure [2] shows the variance of overall capacity R c with 
j3. One can see that R c firstly increases linearly with (3 
and then comes to saturation when (3 is large enough. 
And the saturated value of R c increases with V max . This 
result is in agreement with our common sense that we 
can not improve the overall efficiency only by enhancing 
intersection efficiency, but we should also enhance the 
road condition so that the vehicles can run faster on the 
road. 

In Fig[31 one can see that R c increases with V max until 
a saturation is reached. Therefore, when the intersection 
capacity is large enough, the road condition will be a lim- 
itation for the whole system; while when the road con- 
dition is good enough, the intersection capacity will be 
crucial. Unfortunately, the vehicle speed can not be too 
large in the city. Thus there will be a unavoidable limit 
in improving urban traffic efficiency simply by enhancing 
intersection capacity, given that the network topology is 
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fixed. One should think about other ways, such as adding 
shortcuts, developing subways, employing better naviga- 
tion guidance system for the drivers, and so on. 



IV. SIMULATION OF URBAN TRAFFIC ON A 
SELF-ORGANIZED SCALE-FREE NETWORK 

Recently, works on the centrality of roads in urban 
systems using dual representation [hJ UH [ID, [H, [13, HH 
show that the degree distribution of most planned cities is 
exponential, while it follows a power-law scaling in self- 
organized cities. That is, in most self-organized cities, 
there are some "main roads" with many minor roads in- 
tercrossing with them. 

In this section, we try to simulate the urban traffic on a 
dual graph of scale-free network. To generate the under- 
lying infrastructure, we adopt the well-known Barabasi- 
Albert scale-free network model |8| , in which the power- 
law distributions of degree is in good accordance with 
many real observations. In this model, starting from Too 
nodes fully connected by links with assigned weight wq, 
the system are driven by two mechanics: (1) Growth: 
one node with to links (m < too) is added to the sys- 
tem at each time step; (2) Preferential attachment: the 
probability IF of being connected to the existing node i 
is proportional to the degree fcj of the node 



n, 



(3) 



where j runs over all existing nodes. 

As a remark, here we do not conclude that the BA 
network model exactly describes a self-organized system 
of urban roads. We adopt this model to reflect the fact 
that new roads are usually built to intercross with ex- 
isting main roads. For example, the existing roads are 
often extended to new fields and branch roads are built 
from the extension of these roads. This mechanism can 
lead to the emergence of "main roads" in urban system, 
and it is quite similar to the "growth" and "preferential 
attachment" in BA model. 

Figure U shows the typical evolution of vehicle number 
in the system. The same behavior as in the lattice case 
can be observed. When R < R c , the system is in free 
flow state, and the system will jam when R > R c . 

We first simulate the traffic on a network of N = 100 
nodes (roads) with mo = m = 5. This relatively small 
system can be seen as simulating the backbone of a citys 
urban traffic network. In such a scale-free network, it 
is very important to investigate the effect of navigation 
strategy on the overall capacity. Figure[5]depicts the vari- 
ation of R c with </>. It is shown that R c will be optimized 
at some typical value of </> c . For the case of V max = 100, 
when P is above one, the system's overall capacity will 
be optimized when <f) c — —0.7 and with the R c w 42. 
When p decreases below one, the system's efficiency will 
decrease rapidly. And the optimal value of <fi will increase 
with the decrease of (3. When = 0.3, 4> c = 0.0, implies 
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FIG. 4: (Color online) Typical evolution of vehicle number in 
the scale-free system. The critical generating rate of vehicles 
is R c = 42. 
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FIG. 5: (Color online) The system's capacity vs <j> with 
V m ax = 100. The data are obtained by averaging R c over 
ten network realizations. 



that the best strategy is random-walk. We note that 
4>c < means to use the minor roads first. Therefore, 
if intersection capacity and road condition are the same 
for both main roads and minor roads, the best strategy 
for the whole system is to encourage drivers to use minor 
roads first. 

Figure [5] shows the variation of maximum value of R c 
(the peak value in FigJS]) with the increment of (3 for 
different value of V ma x ■ One can see that R cm ax increase 
first and then come to saturation. And in FigjTl the 
variation of R cm ax with V max is shown. The behaviors 
are similar with the lattice grid case. 

Finally, we try to reproduce the dependencies of aver- 
age velocity and traffic flux on vehicle density, which are 
important criteria for evaluating the transit capacity of 
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FIG. 6: (Color online) The system's maximal capacity Rcmax 
vs f3 with different Vmax- 



70 
60 
50 

X 

ro 

E 40 
o 

or 

30 
20 



BAN=100 
— ■— p=0.3 

— *— p=0.5 
-o-p=1.0 
— a— p=1.5 



50 60 70 80 90 100 110 120 130 140 150 
Vmax 



FIG. 7: (Color online) The system's maximal capacity R c 
vs Vmax with different j3. 



FIG. 8: (Color online) Average velocity vs density. The black 
square line shows the velocity variation when adding vehi- 
cles to the system (increase density), while the red circle line 
shows the velocity variation when drawing out vehicles from 
the system (decrease density) . The density values correspond- 
ing to the sudden transitions are 0.29 and 0.02 respectively. 
The arrows are showing the hysteresis as a guide for the eyes. 
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a traffic system. To simulate the case of constant vehicle 
density, the number of arrived vehicles at each time step 
is recorded, and the same number of vehicles are added 
to the system at the beginning of next step. In FigO the 
velocity-density relation is displayed. Here the vehicle 
density of the system is calculated with p = N c /Ct- The 
velocity firstly decreases gradually with the density. At 
the density of 0.29, the velocity suddenly drops to zero, 
indicating the system enters the jam state. Two branches 
of the fundamental diagram coexist between 0.04 and 
0.29. The upper branch is calculated by adding vehi- 
cles to the system (increasing density), while the lower 
branch is calculated by randomly removing vehicles from 
a jam state and allowing the system to relax after the 
intervention (decrease density). In this way a hysteresis 
loop can be traced (arrows in FigjS]). In the lower branch, 



FIG. 9: (Color online) Flux vs density. 



the velocity keep zero until the density is very low. This 
is because some main roads are congested, thus the ve- 
hicles can not move on these roads, and this state will 
not be alleviated by removing vehicles randomly from 
the system. 

Then we investigate the dependence of traffic flux on 
density. Two kinds of traffic flux are studied: the move- 
ment flux and the arriving flux. The movement flux is 
calculated as the product of average velocity and vehicle 
density. It corresponds to the average number of vehicles 
passing a given spot in the system per time step. Fig- 
ure [HI shows this flux-density relation of the system. The 



6 



10 



- Increase 
Decrease 



? cccmxcaxccamxccccaxn: 



0.0 



0.1 



0.2 0.3 
Density 



0.4 



0.5 



FIG. 10: (Color online) Arriving flux vs density. Arriving flux 
vs density in the cases of increasing and decreasing density. 
Other parameters are the same as in Fig(8] 



arriving flux is calculated as the number of vehicles that 
successfully reach their destination in each time step. In 
Fig llOi the relation of arriving flux vs traffic density is 
shown. In both cases, the hysteresis exist between the 
same values of density as in FigfS] And the maximum 
arriving flux (45 vehicles per step) is corresponding to 
R c = 44 when V max = 100 and (3 = 1.0 as shown in 
Fig© 

The system's sudden drop to the jam state indicates 
a first-order phase transition. The phase transition and 
hysteresis can be explained as follows. According to the 
evolution rules, when a road is full of vehicles, the ve- 
locity will be zero and the vehicles on neighboring nodes 
can not turn to it. So the vehicles may also accumu- 
late on the neighboring nodes and get congested. This 
mechanism can trigger an avalanche across the system 
when the density is high, thus a sudden phase transition 
happen at this point. As for the lower branch, starting 
from an initial congested configuration, the system will 
have some congested roads that are very difficult to dis- 
sipate. These roads will decrease the system efficiency 
by affecting the surrounding roads until all roads are not 
congested, thus we get the lower branch. 



V. CONCLUSIONS AND DISCUSSIONS 

In conclusion, the urban traffic is simulated using a 
model based on dual approach. The model considers 



both the movement of vehicles on the road and the nav- 
igation of drivers in the system. The intersection and 
road conditions are naturally incorporated, and their ef- 
fects on the whole system's efficiency are investigated. 
In a systemic view of overall efficiency, the model repro- 
duces several significant characteristics of urban traffic, 
such as phase transition, hysteresis, velocity-density re- 
lation and flux-density relation. The simulation results 
are a kind of similar to that of the BML model [1] , which 
was proposed in 1992 to simulate the urban traffic on a 
square lattice. But the present model is more general, 
and with more interesting findings. 

As comparing the simulation results on well-planned 
grid and on self-organized scale-free network, one can 
see that the grid network are more efficient than a self- 
organized one. Given that the total capacity (Ct ~ 
57, 600 and 50, 000) is almost the same for the two sys- 
tems considered in this paper, the maximal number of 
vehicles running on the grid is more than that on the 
scale- free network (Fig Q] and HJ, and the overall capac- 
ity R c takes much larger value on the lattice grid (Figf2] 
and [H]). This is in agreement with the previous stud- 
ies showing that homogeneous networks can bear more 
traffic because of the absence of high-betweenness nodes 

When considering the navigation on a scale-free net- 
work, the results are in agreement with previous studies 
that the traffic system will be more efficient by avoiding 
the central nodes [27l ]. 

The work can be extended further in many ways. One 
can modify the model to capture more details in real 
traffic, such as the role of traffic lights, the differences 
between main roads and minor roads. Better navigation 
strategies can also be coined in the dual perspective. The 
resilience of traffic system against road failures will be of 
great importance and research interest. 
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